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ABSTRACT 

We present three new eclipsing white-dwarf / M-dwarf binary systems discovered during a search for tran- 
siting planets around M-dwarfs. Unlike most known eclipsing systems of this type, the optical and infrared 
emission is dominated by the M-dwarf components, and the systems have optical colors and discovery light 
curves consistent with being Jupiter-radius transiting planets around early M-dwarfs. We detail the PTF/M- 
dwarf transiting planet survey, part of the Palomar Transient Factory (PTF). We present a Graphics Processing 
Unit (GPU)-based box-least-squares search for transits that runs approximately 8x faster than similar algo- 
rithms implemented on general purpose systems. For the discovered systems, we decompose low-resolution 
spectra of the systems into white-dwarf and M-dwarf components, and use radial velocity measurements and 
cooling models to estimate masses and radii for the white dwarfs. The systems are compact, with periods 
between 0.35 and 0.45 days and semimajor axes of approximately 2R (0.01 AU). The M-dwarfs have masses 
of approximately 0.35M o , and the white dwarfs are all hydrogen-atmosphere with temperatures of around 
8000K, and have masses of approximately 0.5M Q . We use the Robo-AO laser guide star adaptive optics system 
to tentatively identify one of the objects as a triple system. We also use high-cadence photometry to put an 
upper limit on the white dwarf radius of O.O25R (95% confidence) in one of the systems. Accounting for 
our detection efficiency and geometric factors, we estimate that 0.08%^ ^ (90% confidence) of M-dwarfs 
are in these short-period, post-common-envelope white-dwarf / M-dwarf binaries where the optical light is 
dominated by the M-dwarf. Similar eclipsing binary systems can have arbitrarily small eclipse depths in red 
bands and generate plausible small-planet-transit light curves. As such, these systems are a source of false 
positives for M-dwarf transiting planet searches. We present several ways to rapidly distinguish these binaries 
from transiting planet systems. 



1. INTRODUCTION 

Large numbers of non-eclipsing white-dwarf / main- 
sequence binaries have been disco vered in the Sloan Digital 
Sky Survey and other surv eys (e.g. iReb assa-Manserga s et al.1 
2011; Bian chi et al.1 120071 and references therein). For low- 
mass stars in particular there is a bridge in color between 
white dwarfs and M-dwarfs. The bridge is interpreted as be- 
ing due to rare white-dwarf / M-dw arf binaries, at a rat io with 
respect to single stars of ~ 1:2300 (Smolcic et al. 2004). 

White-dwarf / M-dwarf eclipsing systems are much rarer, 
and almost all have been discovered by searching for white 
dwarfs displaying very deep eclipses of up to several magni- 
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tudes (e.g. iDrake et all 120 101) . These searches find systems 
containing relatively hot (> 12000K) white dwarfs and mid- 
to-late M-dwarfs. The discovery rate o f these systems (e.g. 
IDrake et alj|2l)10t iParsons et alj|201 1 alibi) is increasing with 
the advent of large sky surveys. These binaries survived the 
common-envelope phase of their e volution and many will be- 
come cataclysmic variables (e.g. Neb ot Gomez-Moran et al.l 
1201 11) . and so the properties and number statistics of these 
systems can provide windows into two important areas of 
stellar evolution. Precision measurements of the systems 
allow the determination of the masses and radii of two 
types o f stars for which there are relat i vely few measure- 
ments (iNebot Gomez-Moran et all 120091; iKraus et al.1 1201 It 
iPvrzas et al.11201 1112009b . 

In this paper we present three eclipsing white-dwarf/M- 
dwarf systems discovered during the PTF/M-dwarfs search 
for transiting planets around M-dwarfs. In contrast to most 
known eclipsing systems of this type, the systems detected 
in this survey have optical and infrared emission domi- 
nated by the M-dwarf component and contain relatively low- 
temperature (8000K) white dwarfs and relatively early M- 
dwarfs. The shape of the light curves of the detected systems 
is similar to that expected for transiting giant planets around 
M-dwarfs, in particular in having a flat-bottomed eclipse with 
a depth of 1-20% in red optic al bands. 

The PTF/M-dwarfs survey (lLaw et al.ll201 ll) is a search for 
transiting planets around 100,000 M-dwarfs. The survey is 
perf ormed with the Palomar Transient Factory ( PTF) cam- 
era (iRahmer et alj|2008t lLaw et al.ll2009al l2010al) on the 48- 
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inch Samuel Oschin telescope at Palomar Observat ory, and is 
a Key Project of the P alomar Transient Factory dLaw et al.l 
I2009at iRauet al.ll2009h . The PTF/M-dwarfs survey is de- 
signed to complement o ther M-dwarf transiting planet surveys 
such as MEarth (e.g. ICharbonneau et all 120091: llrwin et al.l 
120101) . the WFCam transit survey (ISipocz et all 1201 lh and 
the M -dwarfs in the Kepler mission target list dBorucki et alj 
I20T1 . by covering a much larger number of M-dwarfs at 
somewhat lower sensitivity. The survey achieves photometric 
precisions of a few percent for ^100,000 targets, and few- 
millimag precision around a subset of »T0,000 M-dwarfs. 
These systems offer much larger transit depths compared to 
solar-type stars, while their very red colors compared to most 
other stars in the field greatly reduce the probability of a 
blended eclipse producing a difficult-to-detect transit false 
positive. 

The three eclipsing systems presented here were originally 
detected as Jupiter-sized planet candidates during the first 
year of operations of the survey. Follow up of the candidates 
showed large color changes during eclipse and very large ra- 
dial velocity signals, suggesting a hidden hot companion. In 
this paper we detail the properties of these eclipsing M-dwarf 
- white-dwarf systems and explore ways to distinguish them 
from true planetary transits. 

The paper is organized as follows: in section|2]we describe 
the PTF/M-dwarfs survey, its precision photometry methods, 
and its target detection strategies, including a new method of 
performing a box-least-squares transit search in parallel on 
GPU hardware. In section[3]we detail the three new eclipsing 
white-dwarf - M-dwarf systems and describe follow up photo- 
metric, low-resolution spectroscopic, and radial velocity ob- 
servations, which are used to determine masses and radii for 
the system components in section|4] Section[5]determines the 
frequency of eclipsing binaries such as these and discusses 
ways to distinguish them from transiting planets. 

2. THE PTF/M-DWARFS SURVEY 

The PTF/M-dwarfs project consists of a transiting planet 
survey on the robotic 48 -inch Samuel Oschin telescope (here- 
after P48), photometric follow-up using the Palomar 60-inch 
telescope, the Byrne Observatory at Sedgwick Reserve (here- 
after BOS), the LCOGT Faulkes-North and Faulkes-South 
telescopes (FTN & FTS), and radial velocity follow-up with 
the HIRES instrument on Keck I. 

The 8-square-degree camera on the P48 telescope allows 
the survey to cover ^3,000 M-dwarfs in each pointing (figure 
Q]), for a total of around 100,000 targets per year at galactic 
latitudes of 20 - 35°. The PTF/M-dwarfs survey typically 
observes several fields with an approximately 20-minute ca- 
dence for 4-5 hours per night. Individual fields are typically 
observed for several months, and observations are performed 
throughout the year. All PTF/M-dwarfs data is acquired with 
a 60-second exposure time in the PTF camera's Mould-R filter 
(similar to the Sloan Digital Sky Survey r filter). Data-taking 
is interleaved with the Palomar Transient Factory supernova 
survey which generally operates on 1-3 days cadences, and so 
scheduling constraints lead to a variety of final cadences for 
PTF/M-dwarfs fields. Table Q] summarizes the specifications 
of the PTF camera and the PTF/M-dwarfs survey. 

The PTF/M-dwarfs survey is sensitive to equal-mass bina- 
ries around all M-dwarfs in the survey fields brighter than 
mp>»:20. In practise, the mass and radius determination and 
follow-up of such faint systems is extremely challenging, so 
we impose a discovery magnitude limit of iiir«18. At that 



TABLE 1 

The specifications of the PTF Camera and the PTF/M-dwarfs survey. 
P48 PTF camera specifications 



Telescope 


Palomar 48-inch (1.2m) Samuel Oschin 


Camera field dimensions 


3.50x2.31 degrees 


Camera field of view 


8.07 square degrees 


Light sensitive area 


7.26 square degrees 


Image quality 


2.0 arcsec FWHM in median seeing 


Filters 


g' & Mould-R; other bands available 


CCD specs 


2KX4K MIT/LL 3-edge butted CCDs 


Plate scale 


1.01 arcsec / pixel 


Readout noise 


< 12 e~ 


Readout speed 


35 seconds, entire 100 MPix array 


PTF/M-dwarfs survey characteristics 


Targets 


Late-K, M and L dwarfs with itir < 1 8 


Survey sky area 


29 square degrees every 2 months 


Target locations 


20° < galactic latitude < 35° 


Targets covered 


~ 12,000 every 2 months 


Observations per night 


5 hours 


Exposure time 


60 seconds 


Cadence 


15-25 minutes 


Observation length 


1-3 months 


Efficiency 


66% open-shutter (slew during readout) 


Saturation Level 


iriR ~ 14 (seeing dependant, in 60s) 


Sensitivity (median) 


iriR«j21 in 60 s, 5cr 




m„»»21.3 in 60 s, 5cr 


Photometric stability 


3 millimag (brightest targets) 




10% (faintest targets) 


Followup 


Photometric: Palomar 60-inch, FTN/S, BOS 




Low-res spectroscopic: Lick Shane-3m 




Radial velocity: Keck I / HIRES 



magnitude, the typical P48 data photometric precision is a 
few percent per datapoint, allowing immediate high-precision 
constraints on the system properties in discovery data. Given 
the saturation and faint cutoff limits of the survey, its effective 
distance ranges are 200-1300 pc for M0 dwarfs, 50-290 pc for 
M5 dwarfs, and 10-70 pc for M9 dwarfs. 

2.1. P48 Data Reduction 

Reduction of the PTF/M-dwarfs data is performed in two 
steps: the standard PTF realtime data reduction software first 
calibrates the images, and then a custom pipeline performs 
source-extraction, association and precision photometry. 

Immediately after observations, PTF data is transferred to 
Lawrence Berkeley National Laboratory where crosstalk cor- 
rections are applied to each chip, standard bias/overscan sub- 
traction is performed, and a superfiat based on recently ac- 
quired data is applied. 

After the calibrated data is transferred to the Dunlap Insti- 
tute for Astronomy and Astrophysics at University of Toronto, 
the PTF/M-dwarfs pipeline extracts sources from the cali- 
brated images, produces an optimal photometric solution, as- 
sociates sources in images taken at different times, and ap- 
plies a range of eclipse-detection algorithms to the resulting 
light curves. We summari ze the syst e m her e; the pipeline is 
described in more detail in lLaw et al.l (1201 lh . 

Initial source extraction is performed on individual CCD 
chip images by SExtractor dBer tin & Ar noutsl 119961) us- 
ing radius-optimized aperture photometry with a locally- 
optimized background. The extracted sources are filtered 
to remove those close to bad pixels, diffraction spikes from 
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Fig. 1. — A PTF Camera image of a survey field centered at 17:28, +57:22 and covering 3.50° X 2.31°. The highlighted points show the 2,851 stars with 
photometrically-estimated spectral types later than K4 and photometric stability better than 5%. The colors of the points correspond to the stellar temperatures, 
with K4 as yellow and late-M-dwarfs as red. The pipeline has removed stars with possible photometric precision problems such as proximity to a bright star or a 
bright ghost image. North is up and East is to the right. 

bright stars, and those that may be affected by nearby sources. 
Heliocentric Julian dates (HJD) are used for time measure- 
ments throughout the pipeline. 

The photometric zeropoints for each epoch are initially es- 
timated based on either SDSS or USNO-B1 photometry for 
bright stars in the field. The pipeline then optimizes the 
zeropoint of each epoch to minimize the median photometric 
variability of all the remaining sources. This first optimiza- 
tion typically improves the long-term photometric stability to 
below the percent level. The pipeline then filters the generated 
light curves, searching for epochs that produce anomalous 
photometry for a large fraction of the sources; those epochs 
are usually those affected by clouds, moonlight or some other 
effect that varies across the images. Typically 0-2% of epochs 
are flagged by this process, and are removed from further 
consideration. A second iteration of variable-source removal 
and zeropoint optimization is then performed. The final ze- 
ropoints are applied to each light curve, along with flags for 
poor conditions, nearby sources that could cause confusion, 
bad pixels, and other problems that could affect the photom- 
etry. Running on a 2.5 GHz quad-core desktop computer the 
pipeline processes a 300-epoch set of 1 1 chips (54 GB of im- 
age data) in less than 24 hours. 

The pipeline typically achieves a photometric stability of 
3-5 millimags over periods of months (figure |2). The photo- 
metric precision is photon-limited for all sources fainter than 
mR =s 16, except in regions of crowding or nebulosity. The 
pipeline has been used for several PTF projects such as the 
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Fig. 2. — The photometric stability achieved in a typical PTF/M-dwarfs field 
in 140 epochs spanning 92 days of observations. Each point corresponds to 
the RMS variability of one of the 23,713 stars in this field with magnitude 
14 < mR < 20. For clarity, 50% of the fainter stars have been removed from 
this plot. The best stability achieved in this field is 4-millimags, with some 
CCD chips and sky regions having somewhat lower stability at 5-7 millimags. 
The few stars with very high variability compared to others at similar mag- 
nitudes are astrophysically varying sources such as eclipsing binaries, RR 
Lyrae stars and other variables. 

open cluster rotation project ( Agii eros et alJl20iTl) . 

2.2. Eclipsing Binary Detection 
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The pipeline produces light curves for 25,000-100,000 stars 
per field. The large number of light curves makes the pro- 
cessing time involved in searching for eclipses an important 
consideration. An initial cut is made on the basis of the 
estimated spectral type of the source, based on photometry 
from the USNO-B1, 2M ASS and (where available) Sloan 
Digit a l Sky Survey dat a dMonet et all 120031: ISkrutskie et alj 
I2006t lYorketdJlioOOl) . The photometric data is fit to up- 
dated versions of the spect ral energy distributions given in 
Krau s & Hillenbrand! |2007), giving an accuracy of approx- 
imately 1 subclass. Sources with late-K, M or L estimated 
spectral types are passed to our eclipse search algorithms. 

2.2.1. High-variability source searches 

Many eclipsing binaries have eclipse signals in the tens- 
of-percent range as well as large eclipse duty cycles. These 
systems have significantly increased photometric variability 
compared to nearby stars of similar magnitude, and so a sim- 
ple variability search can rapidly find them (figure|2]l. We esti- 
mate the value and scatter of the locally expected photometric 
stability as a function of stellar magnitude in 0.1 -mag bins, 
using a sigma-clipped average of the stars detected on each 
chip. Objects that show a more than 2cr increased variability 
compared to the ensemble expectation are flagged for further 
review. Of the systems presented here only PTFEB 11.441 
was detected in this manner; the photometric variations from 
the smaller eclipse depths of the other systems required more 
computationally-intensive algorithms. 

2.2.2. A Parallel Eclipse Search Using Graphics Processing Units 

Systems with smaller eclipse depths and/or longer peri- 
ods and reduced duty cycles may have only slightly in- 
creased photometric variability, necessitating a more sen- 
sitive search . We u se a standard box-least-squares (BLS; 
iKovacs et all (120021) : iTingleyl (120031) ) algorithm to phase the 
light curves at all possible periods and search for a transit-like 
signal. 

This algorithm requires the testing of thousands of peri- 
ods and transit phases and is thus computationally expensive. 
However, the problem is easily parallizable as an arbitrary 
number of light curves and periods can be tested simultane- 
ously. We take advantage of this by implementing the BLS 
algorithm on a Graphics Processing Unit that can perform 
hundreds of computations in parallel. 

Our BLS search is run on an NVIDIA Tesla C2050 Graph- 
ics Processing Unit that contains 448 cores operating at 1.15 
GHz, for a total of 1.03 Tflops in single precision floating 
point arithmetic. The BLS search is coded in CUDA and 
is cal led from the PyCUDA python module (Klockn er et alj 
120091) . 

Perhaps the simplest parallization technique for the BLS 
algorithm is to have each GPU thread perform all the calcu- 
lations for one source; in this way several hundred sources 
could be tested simultaneously. However, the BLS algorithm 
requires at least one access to the full light curve data at each 
test period. The per-thread memory in typical GPUs is too 
small to contain a full light curve, necessitating frequent calls 
to the GPU global memory. These calls are slow, even when 
synchronized across threads, greatly limiting the speed of the 
GPU implementation. 

Instead, we consider sources sequentially, and simultane- 
ously test hundreds of different periods on a single source. 
The processing of each light curve proceeds as follows: firstly, 



the light curve data is read into fast shared memory by each 
thread block (with memory accesses coalesced for speed), 
then each thread picks an untested period, phases the light 
curve at that period and then tests a range of durations and 
phases. 

In testing with typical PTF/M-dwarfs data, our GPU BLS 
algorithm operates approximately eight times faster than a 
multi-threaded program executing on all cores of an Intel 
Core-2 Quad Core CPU running at 2.50 GHz. With the GPU 
approximately 10,000 light curves can be fully searched for 
transit or eclipse events each hour. 

3. DISCOVERIES AND FOLLOW UP OBSERVATIONS 

We discovered three eclipsing M-dwarf / white dwarf bi- 
nary systems in a search of ^45,000 M-dwarfs from the first 
year of PTF/M-dwarf operations. Each system was detected 
as a high-confidence planet candidate, with a 5-20% depth 
flat-bottomed eclipse in R-band with duration consistent with 
a transiting Jupiter-radius planet. Images of the field around 
each target are shown in figure [3] and the P48 detection light 
curves of the systems are shown in figure |4] 

The PTF names of the three detected sources are 
PTF1 J004546.0+415030.0, PTF1 J015256.6+384413.4 and 
PTF1 J015524.7+373 153.8. For brevity we hereafter re- 
fer to the sources by their PTF/M-dwarfs survey inter- 
nal names, which are based on their decimal RA coordi- 
nates: PTFEB1 1.441, PTFEB28.235 and PTFEB28.852 re- 
spectively. 

3.1. Follow up measurements 

After the systems were detected in PTF data as planet candi- 
dates we followed the standard PTF/M-dwarfs follow up strat- 
egy: high-cadence multicolor photometry with the Palomar 
60-inch, BOS, and FTN, along with radial velocity observa- 
tions with the HIRE S spect rograph on Keck- 1 . The multicolor 
photometry (section lS.l.U rapidly revealed that these objects 
had a strongly varying eclipse depth with color (figures [5] 
and|6]l, suggesting we were seeing eclipses between two self- 
luminous objec ts with different temperatures. Low-resolution 
spectra (section [3.1.4t were sufficient to immediately confirm 
a white-dwarf component in one system (PTFEB 1 1.441) but 
required detailed modelling to recover t he w hite dwarf com- 
ponents in the other two systems (section l4~TT i. Radial velocity 
observations were scheduled to allow determinations of the 
white dwarf mass. 

3.1.1. High-cadence, multi-color photometry with BOS 

Multi-color photometry data for each of the targets were 
gathered with the RC Optics 0.8m telescope at Byrne Obser- 
vatory at Sedgwick reserve near Santa Ynez, CA. The BOS 
telescope is equipped with a Santa Barbara Instrument Group 
STL-6303E camera utilizing a 3k x 2k Kodak Enhanced KAF- 
6306E CCD, with a 14.7' x 9.8' field of view and a pixel scale 
of 0.572" per pixel (2x2 binning). We observed individual 
eclipses in a single color with SDSS i', SDSS g', or Astrodon 
Photometries UV-blocked clear filters. 

The images were reduced using standard routines for bias 
subtraction, dark current subtraction, and flat-field correction. 
We extracted fluxes for all stars in the frame using aperture 
photometry routines in PyRAF. Relative flux light curves were 
produced by dividing the flux of the target star by the summed 
flux from several comparison stars in each image. Julian 
dates of mid-exposure were recorded during the observations, 
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Fig. 3.— Images of the newly-discovered systems taken with the PTF camera (left-to-right: PTFEB 11.441, PTFEB28.235 and PTFEB28.852). North is up and 
East is to the left; each image shows a 400x400 arcsecond cutout from the 3.50 X 2.31 degree PTF camera field. 
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Fig. 5. — Faulkes-North (z-filter) and BOS (g-filter) follow-up photome- 
try of PTFEB28.852. The eclipse is detected at high significance in both 
bands, but is approximately lOx deeper in the g filter, suggesting the sec- 
ondary eclipse of a hot body. 
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Fig. 4. — The discovery P48 light curves for the three eclipsing systems (top- 
to-bottom: PTFEB1 1.441, PTFEB28.235, PTFEB28.852). The colouration 
of the points corresponds to the time the datapoint was measured; black are 
the oldest points and bright yellow are the newest. PTFEB 11.441 was ob- 
served during a high cadence run targeted at M3 1 , hence the large density of 
points all taken on a single night. 

and later converted to ba rycentric Jul i an dat es using the on- 
line tools described by lEastman et al.l (120101) . Aperture sizes 
were optimized by minimizing the scatter of the resulting light 
curves, and ranged from 3" to 6". Exposure times were typi- 
cally 300s. 

Several eclipses of PTF1 1.441 were observed in short- 
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Fig. 6. — Dunlap Institute Arctic Telescope photometry of PTFEB 11.441, 
showing the varying eclipse depth with wavelength. The g, r and i filters 
(with midpoints shown in the figure) were alternated in sequence to produce 
a near-simultaneous multi-color light curve. To guide the eye, each color has 
a fit to the light curve using the system's eclipse parameters derived in section 
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cadence mode at BOS. We used an Astrodon Photometries 
UV-blocked clear filter in order to maximize the signal to 
noise ratio (SNR) per unit time. We used 3x3 binning (0.858" 
/ pixel), and set the camera to only readout a small subsection 
of the CCD encompassing the target star, and three nearby 
stars of similar brightness to be used as comparison stars. This 
reduced the readout plus dead time of the instrument from 
« 10s to a;5s, and allowed us to achieve a photometric noise 
rate (PNRHof ~ 1.6%/min. 

3.1.2. Faulkes-North Photometry 

PTF28.852 was observed photometrically on 2010-11-10 
and 2010-11-16 using the 2.0 m Faulkes Telescope North 
(FTN, Maui, Hawai' i) operated by LCOGT. In both cases, the 
Spectral CCD imageiFH was used along with a Pan-STARRS- 
z filter. The Spectral instrument contains a Fairchild CCD486 
back-illuminated 4096 x 4096 pixel CCD which was binned 
2x2 giving 0.303" pixels and a field of view of 10' x 10'. Ex- 
posure times were 130 & 150 seconds respectively on the two 
nights. 

The frames were pre-processed using standard techniques 
for bias subtraction and flat-fielding; dark current subtraction 
was not performed as it is negligibly small for this instrument. 
Object detection and aperture photometry were performed us- 
ing the DAOPHOT photometry package within the IRAF en- 
vironment. The aperture sizes used were 7 and 5 pixels in 
radius on the two nights. Differential photometry was per- 
formed relative to 5-6 comparison stars within the field of 
view. 

3.1.3. Dunlap Institute Arctic Telescope Photometry 

The Dunlap Institute Arctic Telescope is a 20-inch robotic 
telescope currently undergoing testing at the New Mexico 
Skies observatory at Cloudcroft, NM. Once testing is com- 
plete, the telescope will be based at the PEARL research sta- 
tion on Ellesmere Island at a latitude of 80 degrees North, 
where it will perform a search for transiting planets. The tele- 
scope is equipped with a 16-megapixel Apogee U16M camera 
with a 34' x 34' field of view. 

We observed PTFEB 11.441 on the night of October 23 
201 1 . 120-second exposure g, r, and i-band images were taken 
in sequence throughout an eclipse of the target. After standard 
calibrations, differential pho tome try was performed using the 
pipeline described in section |2~T1 

3.1.4. Low-resolution Spectra 

We obtained low-resolution optical spectra of PT- 
FEB28.235 and PTFEB28.852 on 2010 Nov 30 UT, and PT- 
FEB 11. 441 on 2011 Oct 25 UT Spectra were acquire d with 
the dual-arm Kast spectrograph dMiller & Stond [T993) on the 
3-m Shane telescope at Lick Observatory. The spectra used 
a 2 arcsec wide slit, a 600/4310 grism on the blue side and a 
300/7500 grating on the red side, yielding FWHM resolutions 
of * 4 A and w 10 A. All observations were aligned along the 
parallactic angle to reduce differential light losses. 

All spectra we re reduced using standard techniques (e.g., 
iFolev et aTll2003l) . Routine CCD processing and spectrum ex- 
traction were completed within IRAF, and the data were ex- 

10 Photometric noise rate, calculated as RMS / V(r), where RMS is 
the scatter of the out-of-eclipse section of the light curve and T is the me- 
dian number of cyc les (exposure time and dead time) per minute. See 
Shporer et al. 1 2010b) for a more detailed description. 

1 1 http://www.specinst.com 
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tracted with the optical algorithm of iHornd (119861) . We ob- 
tained the wavelength scale from low-order polynomial fits to 
calibration-lamp spectra. Small wavelength shifts were then 
applied to the data after cross-correlating a template sky to the 
night-sky lines that were extracted with the supernova. Using 
our own IDL routines, we fit spectrophotometric standard-star 
spectra to the data to flux calibrate our spectra and to rem ove 
telluric lines (IWade & Horndll988tlMatheson et aT]|2000l) . 

3.1.5. Keck/HIRES Radial Velocity Observations 

We obtained high-dispersion spectra for the whte-dwarf / 
M-dwarf systems using the High-Resolution Echelle Spec- 
trometer (HIRES) on the Keck-I 10m telescope. HIRES is a 
single-slit echelle spectrograph permanently mounted on the 
Nasmyth platform. All observations were performed using 
the red channel, and span a wavelength range of 4300-8600 
angstroms. All were obtained using the C2 decker, which 
yields a spectral resolution of R ~ 45000. We processed our 
HIRES data using the standard pipeline MAKEE, which auto- 
matically extracts, flat-fields, and wavelength-calibrates spec- 
tra taken in most standard HIRES configurations. In tabled 
we list the epochs and exposure times for all of our HIRES 
observations. 

In order to remove small drifts in the wavelength calibra- 
tion over the course of the night, we cross-correlated each 
spectrum's telluric features at 7600 angstrom s with the tel- 
luric f eatures for standard stars selected from iNidever et all 
(2002), placing all observations into a common frame set 
by the Earth's a tmosphere. As we have previously shown 
dKraus et alJl201 ll) . this reduces the systematic uncertainty in 
the measured velocities to ~0.3 km/s. We then measured the 
position and equivalent width (EW) of the H-alpha line at each 
epoch. The results are summarized in table[2] 

3.1.6. Robo-AO Laser Guide Star Adaptive Optics Observations 

Robo-AO is a visible and near-infrared laser guide star 
adaptive op tics system specifica lly engineered for 1-3 m class 
telescopes (Bar anec et alJl20lH) . The Robo-AO system com- 
prises an ultraviolet Rayleigh laser guide star, an integrated 
adaptive optics and science camera system, and a robotic 
control system. The system currently incorporates both an 
electron-multiplying CCD (EMCCD) and an InGaAs infrared 
array camera for imaging. 

Robo-AO obtained images of PTFEB 11.441 and PT- 
FEB28.235 during a Robo-AO commissioning run at the 
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Fig. 7. — Robo-AO laser guide star adaptive optics images of (left-to-right) PTFEB1 1.441, PTFEB28.235, and an m, ~ 13 binary taken in the same part of the 
night to provide a comparison point spread function. All images are in a 600nm long-pass filter, are 4 arcseconds square, have the same rotation, and are displayed 
with linear scaling. To improve the signal-to-noise for these faint targets, the PTFEB1 1.441 and PTFEB28.235 images have been smoothed by a Gaussian filter 
with a FWHM corresponding to the diffraction-limited resolution of the telescope. The apparent vertical smearing of the PTFEB 11.441 image is due to the 
shift-and-add algorithm occasionally locking on to the secondary component of the binary (Law 2007). 



Palomar 60-inch telescope on the night of 2011 Nov 18 UT. 
Robo-AO was operated without tip-tilt correction, instead re- 
lying on post-facto shift-and-add processing of the individual 
frames. We used a long-pass filter with a 600nm cut-on to 
obtain increased signal compared to a bandpass filter, along 
with relatively long 1 -second exposure times to compensate 
for the faintness of the targets (hir > 16). We used the Lucky 
Imagin g plus adaptive optics pipeline described in lLaw et al.1 
(2009b) to perform the image alignment, along with frame se- 
lection at the 10% level to boost the imaging resolution. After 
processing, images with approximately 0.2" full-width half- 
maximum (FWHM) point spread functions were obtained on 
these targets during ~1" seeing conditions. These observa- 
tions were limited by residual tip-tilt errors during the 1- 
second exposures; similar bright targets observed the same 
night with 33-ms exposures had a diffraction-limited angular 
resolution of «0.13 arcsec FWHM (figure[7]i. 

The Robo-AO imaging reveals a bright companion to PT- 
FEB1 1.441, located at 0.43+0.06 arcsec and at a position an- 
gle of 165 ±3 degrees East of North. 

3.2. The nature of the detected objects 

Multi-color through-eclipse photometry of the three de- 
tected systems showed that the eclipse depths were much 
larger in blue bands than in red bands. The sharp ingress 
and egress of the eclipses suggest a small body, while the 
flat-bottomed shape argues against that body eclipsing the M- 
dwarf, as limb-darkening would then be apparent. Radial ve- 
locity observations showed a very large amplitude of more 
than lOOkm/sec. These observations suggest the presense of 
a small, hot, massive object undergoing total eclipses in these 
systems. The white-dwarf / M-dwarf binary scenario is con- 
sistent with all these observations. 

The properties of the systems are summarized in table [3] 
We detail each system below. 

3.2.1. PTFEB1 1.441 

This system shows a possible companion in Robo-AO im- 
ages. If it is associated, the companion is quite wide, at ap- 
proximately 80 AU separation. There are no other sources de- 
tected in the Robo-AO image, which covers 1600arcsec 2 . On 
this basis, the probability of an unassociated star appearing 
within 2 arcsec of our target is less than 1 %. Definitive confir- 
mation will require colors and common-proper-motion mea- 



surements, but for the purposes of the following sections we 
tentatively associate the companion star with PTFEB 1 1.441. 

The two wide components of the system have roughly equal 
brightness and the low-resolution spectrum of the system is 
well fit by a combination of a white-dwarf an d an M-dwarf 
spectrum without further components (section 14. Il l, suggest- 
ing that spectral types of the two M-dwarfs are very similar. 
Wide M-dwarf binaries make up a small fraction o f the binary 
popul ation, but many examples have been found (Dhit al et alJ 
2010), and a very large fraction of the components of the 
widest binaries have clos e companions with roughly equal 
masses (lLaw et alJ l2010bl) . However, we note that the M- 
dwarf near the white dwarf in PTFEB 11.441 survived the 
common-envelope phase of the system's evolution, and so it is 
unlikely that this system formed in the same manner as other 
known wide M-dwarf binaries. 

PTFEB 11.441 lies in a PTF field targeted at M31 and 
is located approximately one degree from the center of the 
galaxy. Of the three systems, PTFEB1 1.441 shows the great- 
est amount of emission from the white dwarf, with the M- 
dwarf being overwhelmed at wavelengths < 550nm. The 
light curve shows only very-low-level out-of-eclipse variabil- 
ity. The target is detected by GALEX with a near-UV magni- 
tude of niNuv = 17.35, with no far-UV detection. 

3.2.2. PTFEB28.235 

The P48 light curve of this system shows out-of-eclipse 
variability at the 7-10% level, almost as strong as the eclipse 
itself in R-band. The variability phasing (maximal when the 
white dwarf is in front of the M-dwarf) suggests that the vari- 
ability is primarily an irradiation effect. The 5% eclipse depth 
is the smallest R-band depth of the three systems. The system 
has a GALEX Near-UV detection at m NUV = 21 .6 ± 0.4 (with 
no far-UV coverage). 

3.2.3. PTFEB28.852 

This system is very similar to PTFEB28.235, with similar 
variability levels and a similar light curve shape. The system 
has a GALEX far-UV detection at m FUV = 21.0 + 0.1 (with 
no near-UV coverage). 

4. MASS AND RADIUS MODELS AND MEASUREMENTS 

The optical light in these systems is dominated by the M- 
dwarfs, and the large radius ratios between the components 
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TABLE 3 

Measured and Estimated System Properties 



Name RA Dec Period To SpTmd Phot.SpTMc T^d Dist./pc Mmd Rmd 

(HJD) (±0.5) (±0.5) (±500K) 

PTFEB 11.441 00 45 46.0 +41 50 30.0 0.35871d±0.00005d 2455438.3165 M3 M3~5 8500K 180±45* 0.35±0.05 M 0.33±0.05 R 

PTFEB28.235 01 52 56.6 +3844 13.4 0.3861 ld±0.0001 Id 2455460.1116 M3 M3 8000K 200±5() 0.35±0.05 M 0.33±0.05 R s 

PTFEB28.852 01 55 24.7 +37 31 53.8 0.46152d±0.00009d 2455530.7335 M3 M2 8500K 260±7() 0.35±0.05 M 0.33±0.05 Ro 

Note. — The properties of the detected systems, based on photometry and low-resolution spectroscopy. SpT^D is the M-dwarf spectral type derived from 
low-resolution spectra, and Phot.SpT M D is the spectral type estimated from infrared colors. The white dwarf properties are estimated as described in the text. The 
distance estimate for PTFEB 11. 441 is starred because it is likely to be an underestimate because of its possible companion, which is unresolved in the 2MASS 
photometry used to generate the distance estimate. 



produce very small (millimag-level) primary eclipses. For 
these reasons, we were able to measure precise eclipse shapes 
only for the white-dwarf occultation events, and radial ve- 
locities for the M-dwarf components. However, the low- 
resolution spectra allow us to estimate the mass and radius 
of the M-dwarf component of the system and thus estimate 
the properties of the white dwarf component. Finally, high- 
cadence photometry of the ingress and egress of the white- 
dwarf occultation events allowed us to set stringent upper lim- 
its on the radius of the white dwarf in PTFEB 1 1.441. 

4.1. Primary and secondary spectra 

We model the low-resolution spectra of our targets as a 
combination of a M-dwarf and a white dwarf spectrum. Tem- 
plate M-dwarf spectr a are taken from the HILIB stellar flux 
library (Pickles 1998) and have steps of a single spectral sub- 
class. The white dwarf model£j were kindly provided by 
Detlev Koester and are on a grid with 25 OK steps between 
6000K and 10000K (and longer steps up to 20000K), and 0.25 
dex steps between log g = 6.0 and log g = 9.5. 

We simultaneously fit the white dwarf and M-dwarf spec- 
tra using a down hill simplex algorithm and a bootstrap (e.g. 
iPress et aT] 1 19921) method to find the uncertainties in the fit. 
Figure [8] shows the fitted models and each component of the 
binaries with the model for the other component subtracted. 
In each case subtraction of the template M-dwarf spectrum 
revealed a (high-noise) DA white dwarf spectrum with clear 
Balmer lines. 

PTFEB 1 1 .441 is well-fit by a simple combined white-dwarf 
and M-dwarf spectrum, suggesting that the third component 
in the system is also an M-dwarf, and that the two M-dwarfs 
have similar spectral types. In the following sections, we use 
the measured spectral type for the M-dwarfs in PTFEB 1 1 .44 1 , 
assuming that it is the spectral type of the M-dwarf closest to 
the white dwarf. 

4.2. Reddening 

The M-dwarf spectral types measured from the two- 
component fit to the system spectra closely agree with those 
estimated from infrared photometry of the systems (table [3). 
This suggests that reddening is not a significant problem for 
the M-dwarf spectral types and thus mass and radius estima- 
tions. 

However, we note that reddening will affect the estimated 
white dwarf temperature more severely as its flux is mostly in 
the blue end of the spectrum. The estimated galactic extinc- 
tion along the line of sight to PTFEB1 1.441 is E(B-V)=0.02 
dSchlegel et al.l lT998). and will be lower for the target itself 
which is located at a distance of « 400pc. The other two 

12 IKoe ster (2010). Balmer lines in the models were calculated with the 
modified Stark broadening profiles of Tremblav & Bergeron (2009), kindly 
made available by the authors. 



targets are located in regions of higher line-of-sight extinc- 
tion, at E(B-V)=1.5 and 1.1 mag for PTFEB28.235 and PT- 
FEB28.852 respectively. The actual extinction for the two 
targets is, however, likely to be much lower as they are lo- 
cated at distances of only ^200pc. This is supported by sev- 
eral different arguments: the very good match of the M-dwarf 
model spectrum to the observed spectrum; the match between 
the infrared-color-estimated spectral types and the spectrally- 
measured spectral types; and the GALEX detections of both 
objects, which require a small near and far-UV extinction. 

4.3. M-dwarf Properties 

We fit the M-dwarf masses, radii and distances in the fol- 
lowing manner: we first fit sur vey photometry of the system 
with the stellar SEDs listed in iKraus & Hillenbrand! d2007l) . 
yielding a photometric spectral type and distance. Despite 
the availability of other photometric data such as USNO-B1, 
we restrict our fits to 2MASS J, H and K colors that are 
unlikely to be affected by the light fr om the white dwarf. 
We assume a 25% distan ce uncertainty ( Kra us & Hillenbrand! 
120071; lLaw et ai]l2010al), and ±1-2 spectral subclass uncer- 
tainty !&MiI&!Hinenbran3|2003); the SED estimates for 
the spectral types are consistent with those derived from the 
two-component fits to our low-resolution spectroscopy. We 
also attempted to determine the M-dwarf sp ectral type from 
the Ti05 narrowband spectroscopic index (Gizis [19971) . but 
found implausibly high values of the index, suggesting possi- 
ble contamination by white dwarf emission. In the following 
sections we adopt the spectral types estimated from the two- 
component fits to the low-resolution spectra, as that method 
yields the lowest spectral type uncertainty. 

We estimate the M-dwarf masses from th e spect ral type vs. 
mass calibrations detailed in Delfo sse et alj (|2000), assuming 
they are on the main-sequence and have solar metallicity and 
age. All three of our targets have a measured spectral type 
of M3, giving an estimated mass of 0.35±0.0 5 M n . We note 
that a s mall fraction of the systems measured in Delfos se"et alj 
(2000) have much lower masses at the M3 spectral type, for 
reasons that are still unclear. To estimate the masses in an 
alternate manner, we estimated the M-dwarf effectiv e temper- 
atures from the relation described in Luhman ( 1999). We then 
combined those esti mates with the 5GY r isochrones of solar- 
metallicity stars in Baraf fe"et alj 0998) to estimate the stellar 
masses. These relations also predict a much lower mass for 
our M-dwarf targets, of O.23g l O 4 6 M . In the following sec- 
tions w e adopt the higher masses measured in Delfosse et al. 
(2000). We, however, note that two-component RV mea- 
surements for our systems giving masses for each component 
would help understand the mass / spectral type relation in this 
mass range. 

We estimate the M-dwarfs' radii from a fit to the eclipsing- 
binary-derived mass/radius relation. We note that M-dwarfs 
in close binaries appear to be biased towards larger radii com- 
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Fig. 8. — Spectra for the three targets: in order from top to bottom, PTFEB1 1.441, PTFEB28.235, PTFEB28.852. Left: the observed spectrum (red); the 
M-dwarf (green) and white-dwarf (blue) spectra simultaneously fit to the spectrum, and the fitted combined spectrum (black). Right: separated spectra - the 
observed spectrum with each model spectrum removed in turn. The Balmer lines in the white dwarf spectra are labelled. 



pared to systems with wider orbital radii (iKraus et al.112.01 ll) . 
We adopt a mass/radius relation for close (< 1 da y period) bi- 
naries based on the measurements summarized in Kra us et all 
d20Tlh : 



TABLE 4 
Spectroscopic orbit properties 



R = 1.096M - 0.052 for 0.1 M Q < M < 0.7M f 



(1) 



Name 

PTFEB1 1.441 
PTFEB28.235 
PTFEB28.852 



Ki (km/sec) 

-166±15 
-186±17 
154±5 



Vo (km/sec) 

-29±26 
16±15 
13±3 



Mass function (F) 

0.14±0.04 
0.26±0.07 
0.17±0.02 



where M is the mass in solar masses and R is the radius in 
solar radii. The radii in this relation are approximately 5-10% 
larger than both those measured for lar ger orbital rad i us sys - 
tems and theoretic al expectations (e.g. Bar affe et al.l (1998); 
IKraus et ail (1201 1|) ~). The results are summarized in table [3] 



4.4. White Dwarf Properties 

We measure the white dwarf masses from the measured ra- 
dial velocity curves, assuming the M-dwarf masses derived 
above. We fit a simple model to the velocity curves (fig- 
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TABLE 5 
White dwarf properties 



Name M WD K WDmode i J*WD.spec. K WD,direct 

PTFEB 11.441 0.51±0.09M o 0.014±0.002R o 0.010R o <0.025R o (95%) 
PTFEB28.235 0.65±0.1 1M 0.012±0.002R o 
PTFEB28.852 ().52±0.05M G 0.014±0.001R o 

Note. — Derived masses and radii for the white dwarfs in the eclipsing 
binary systems. Mwd is the stellar mass; RwD.model is the radius derived 
from the mass and the Bergeron cooling models; RwD,spec. is the radius 
derived from direct fits of the Koester (2010) model spectra Balmer lines 
to the white dwarf spectrum; RwD.direct i s me 95%-confidence upper limit 
measured from the eclipse ingress and egress shapes. 

ure [9] and table 0), assuming zero eccentricity because of 
the short circulari zation timescales f or these compact systems 
(for example, see lDevor et"aTll2008l) . PTFEB28.235 and PT- 
FEB28.852 have sufficient numbers of RV points to allow an 
eccentric orbit fit, but the fitted eccentricities were consistent 
with being zero. The white dwarf masses (Mwd) are related 
to the fit mass function (F), the M-dwarf masses (Mmd) and 
the inclination (i) by: 



F = 



(2) 



Fig. 9. — Radial velocity orbit fits for each of the targets. Each binary shows an anomalous measurement around the quadrature points, which may be due to 
the white-dwarf's Ha absorption line interfering with the measurement of the M-dwarf Ha position. A simple radial velocity orbit is consistent (reduced^ 2 a* 1) 
with all the other points, but inconsistent with the anomalous points, and so we remove those points to avoid unduly biasing the fits. 

We fit the light curves with a simple two-spherical-body 
model. The flat-bottomed light curves show that the white 
dwarf is completely ecl i psed, allowing a simple model based 
on the iMandel & Agoll J2002) transiting-planet model to be 
used. We modified the model to allow both stars to be moving 
relative to a common center of mass, to allow both bodies to 
be self-luminous, and to include a linear limb-darkening term 
for the white dwarf. 

Inspection of the light curves shows that the white dwarf 
ingress and egress is at best marginally resolved, and so we 
only attempt to estimate an upper limit for the white dwarf 
radius. We obtain a maximal estimate of the white dwarf ra- 
dius when the inclination is assumed to be 90°, and so we fix 
the inclination to that value in our modelling. For the white 
dwarf, we conservatively use a linear limb-darkening coeffi- 
cient of 0.5; we tested a range of coefficients and found that 
the limb darkening has no significant effect on the radius es- 
timates, because the eclipse ingress and egress are barely re- 
solved. We fix the M-dwarf radius and the stellar mass ratio 
q during the fits, although we marginalize over the estimated 
1 -sigma ranges of both parameters. 

We vary the semimajor axis of the system, the white dwarf 
radius, and the flux from both components (and hence the 
contribution from the system's unresolved wide companion). 
When the inclination is fixed at 90 degrees, the best fit semi- 
major axis value is within the errorbars of our measurements 
based on the eclipse timing and the estimated M-dwarf radius, 
although those estimates suggest a somewhat lower inclina- 
tion may be preferable (which would lead to a smaller white 
dwarf radius). We estimate the range of fitted parameters and 
a 95% upper limit on the white dwarf radius using the Boot- 
strap algorithm. 

This procedure puts the upper-limit white-dwarf radius in 
the PTFEB 11.441 system as <O.O25R at 95% confidence. 
This value admits the radii estimated from cooling models 
and log g fits to the white-dwarf's Blamer lines. The fitted 
semimajor axes are around 10% larger than those expected 
from the period of the system and the estimated masses of the 
components, suggesting we have slightly overestimated the 
mass of at least one component, most likely the M-dwarf in 
the system. 

4.6. M-dwarf Activity 

All three systems show evidence of active M-dwarfs, in- 
cluding strong Ha emission and variability on a variety of 
timescales. Non e are detected in t he ROSAT All-Sky Faint 
Source Catalog (IVoges et al.1 120001) . although the M-dwarfs 
would have to be exceptionally active to be detected at these 
distances (e.g. lLaw et al.| [2008). 



M WD sin i 
(Mmd + Mwd) 2 



which can be solved as a cubic equation. The results are sum- 
marized in table [5] The inclination has only a weak effect on 
the derived white dwarf masses; we marginalize it across all 
ranges allowed by the eclipse geometry. 

We estimate the white dwarf radii in two ways: firstly, mod- 
els of the white dwarf effective temperature, mass and surface 
gravity cooling curves, and secondly a direct fit to the white 
dwarf's eclipse light curve. 

Our model radius estimates are base d on the Bergeron 
cooling curves for DA whi t e dwa rfs (Holberg & Bergeron 
2006t iKowalski & Saumonl l2006t iTremblav et alJ 1201 It 
Berge ron et al.l [20 1 1). for which we derive the allowed sur- 
face gravity range from the white dwarf mass range measured 
above. The effective temperature has only a small effect on 
the surface gravity estimates and so we fix it at the best-fit 
value from the low-resolution spectra. 

4.5. Light curve models 

We obtain a direct upper limit for the white dwarf radius in 
the PTFEB 1 1.441 system from high-cadence eclipse data ac- 
quired with the BOS telescope (figure [T0l. Using a clear filter 
to boost the signal levels, we were able to obtain 10-second- 
cadence data with »5% photometric stability, sufficient time 
resolution to give stringent constraints on the ingress and 
egress timing of the system. 



Three New Eclipsing White-dwarf - M-dwarf Binaries 



11 



1.15 




0.75 



0.00 



0.50 



1.00 
Time / hours 



1.50 



2.00 



Fig. 10. — 10-second-cadence observations of PTFEB11.441 from the BOS 
telescope (black points). A clear filter was used to increase the signal from 
the target. The best-fit eclipse model is shown in red. 
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Fig. 11. — A flare event during a high-cadence observation of PTFEB1 1.441. 
The green line shows a 30-second moving average of the 10-second-cadence 
photometry (black points). The red line shows the best fit eclipse light curve. 

The positioning of PTFEB 11.441 near M31 allows more 
archival data to be found, and we identify the system with 
the Chandra X-ray sourcf] CXO J004545.8+415029 (2.8 x 
10 _13 erg/ s/cm 2 ) located within lcr of the PTF position for the 
system. 

We also identified a flare event during one high-cadence ob- 
servation of PTFEB 1 1.441 (figure ITTb. The observation was 
taken without a filter to boost our signal-to-noise, limiting the 
detailed interpretation of the event. The shape of the event is 
as expected for an M-dwarf flare, with a rapid raise followed 
by an approximately exponential decay. 

4.7. Out-of-eclipse variations 

The light curves show out-of-eclipse sinusoidal variabil- 
ity correlated with the orbit. The overall shape is consistent 
with irradiation of the side of the M-dwarf facing the white 
dwarf, where the flux maximum is expected half an orbital 
phase from the occultation. Other mechanisms that may af- 
fect the out-of-eclipse light curves of such compact systems 
are tidal ellipsoidal distortion and the beaming effect (e.g., 
IZucker et al.l200 7; Shpor eret alj2010al) . but simple estimates 
show that both effects are expected to be at the 1 % level or be- 
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Fig. 12. — The cumulative number of M-dwarf targets searched by the 
PTF/M-dwarfs survey in the first year of operations, as a function of achieved 
photometric stability. M-dwarfs in the PTF/M-dwarf fields are selected and 
confirmed on the basis of 2MASS, USNO-B1 and (where available) SDSS 
colors, along with proper motions determined from 2MASS and USNO-B1 
positions. We require a high-confidence photometric color fit, along with a 
2<x proper motion detection to rule out giant interlopers. 

low for the binaries presented here and so can not be identified 
in the available data. 

PTFEB 11.441 shows much lower levels of out-of-eclipse 
variations than the other two systems. This is most likely due 
to light from its possible companion diluting the variations 
(and eclipse depths) in unresolved photometry of the system. 

5. DISCUSSION AND CONCLUSIONS 

The discovery of only three white-dwarf / M-dwarf bina- 
ries in the PTF/M-dwarfs survey suggests a low incidence for 
such systems. This phase of the survey has covered «45,000 
M-dwarfs with sufficient precision to detect short-period sys- 
tems with similar eclipse depths to those shown here (fig- 
ure fT2l . The survey's detection efficiency for systems with 
fractional-day periods is near 100% (figure \13[ . The geo- 
metric probability of eclipse is 5-15% for an early-M-dwarf 
primary in a roughly equal-mass half-day-period binary with 
a white dwarf. Taken together, the detection of three systems 
implies that 0.08%!° ^* (90% confidence) of M-dwarfs are in 
short-period post-common-envelope white-dwarf / M-dwarf 
binaries where the optical light is dominated by the M-dwarf. 

The periods of the detected systems are all longer than most 
known eclipsing M-dwarf / white-dwarf systems (figure [T4"ii. 
Our survey is designed to detect longer-period transiting exo- 
planets, so this is not in itself surprising. However, our lack of 
detections at shorter periods, despite our near- 100% detection 
efficiency for such systems, suggests that binaries including 
these relatively low-temperature white dwarfs are preferen- 
tially found at relatively large orbital radii. 

These eclipsing binaries appear as M-dwarfs when selected 
by optical and infrared colors, and they display transit-like 
light curves which can have arbitrarily small depths. They 
appear to be at least three times more common than transit- 
ing giant planets in the PTF/M-dwarfs survey. These proper- 
ties make these systems an important false-positive for current 
and future transiting planet searches around M-dwarfs. Fortu- 
nately, they can be distinguished from true transiting planets 
with small amounts of extra data, using the methods detailed 
below. The methods are ordered by difficulty, starting with 
those which require only the discovery light curves. 



13 from the Chandra Source Catalog (Evans et al. 2010) 
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Fig. 13. — The eclipse detection efficiency for a typical PTF/M-dwarfs ob- 
serving scenario, based on simulations including the actual observing ca- 
dences of the survey. The detection efficiency is the probability that at least 
two eclipses will be observed from an eclipsing target with a specified period, 
after observing window and weather effects have been taken into account. 
The low efficiencies around integer numbers of days at short periods are due 
to the likelihood of an eclipse always occurring in daytime; the same effect al- 
lows some long period binaries to be picked up if their eclipses always occur 
at night. 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



CD 

cr 
CD 







t- in 


C\J 


co 


LO 


eg 


CO 


oci 


CO 


t- CM 


C\J 


CO CO 


CO 


LU LU 


LU 


LL LL 


LL 


l-l- 


1- 


D_ Q_ 


Q. 


•• 


• 



0.1 0.2 0.3 0.4 
Period /d 



0.5 



0.6 
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IDrake et all (20T0t) , which were all found in a single survey and so offer a 
useful comparison sample. 

compared to a M-dwarf/planet system leads to a shorter 
eclipse; high-SNR and high-cadence light curves are 
however required to distinguish this case from a lower- 
inclination system. 

2. Eclipse shape. The eclipse of the white dwarf by the 
M-dwarf provides an exactly-fiat-bottomed eclipse out- 
side the short ingress and egress periods. With high- 
precision photometry this can be distinguished from 
a true transiting planet which will show the effects of 
limb darkening as it passes in front of the M-dwarf. 

3. Out-of-eclipse variability. All our systems show out-of- 
eclipse variability at levels which are orders of magni- 
tude greater than t hat expected for planetary systems 
(see iDrakd d2003l) for details of a similar selection 
method for planetary transits around solar-type stars). 

4. UV emission. Two of our new systems showed UV 
emission and the third has possible emission. Where 



data is available, the presence of UV (or even u-band) 
emission suggests the existence of something other than 
an M-dwarf in the system. 

5. Low-resolution spectroscopy. Like the three targets 
presented in this paper, a low-resolution spectrum could 
be decomposed into white dwarf and M-dwarf compo- 
nents. Faint, low-temperature white dwarfs may be, 
however, hard to detect with this technique. 

6. Multi-color photometry. These systems show a strong 
variation in eclipse depth with wavelength (depending 
on the relative temperatures of the white dwarf and M- 
dwarf and their ratio of radii). 

7. Radial velocities. These systems have radial veloc- 
ity amplitudes at least 500x larger than expected for a 
planet with the same period. Just two RV few-km/sec- 
precision datapoints are sufficient to discriminate these 
systems from transiting planet systems. 

Of these methods, mulitcolor photometry through eclipse is 
probably the most time-efficient method of detecting systems 
like this. As large scale transiting planet surveys of M-dwarfs 
start up and continue, many more systems in this interesting 
parameter space are likely to be found. 

Follow-up precision photometry and radial velocities will 
allow direct measurements of the masses and radii of all 
components of these systems. As the white dwarf tran- 
sits across the disk of the M-dwarf the transit depth is ex- 
pected to be around 1 millimag, although lensing by the 
white dwarf (e.g., [Marsh 2001) will make the transit shal- 
lower than expected from only geometric considerations (e.g., 
ISteinfadt et al.l l2010). High-cadence and high-precision pho- 
tometry of the white dwarf eclipse (occultation) ingress and 
egress may be the best approach to directly measure the 
white dwarf radii. Furthermore, small asymmetries in the 
ingress and egress li ght curve due to the photometric R ossiter- 
McLaughlin effect (Shp orer et al.ll2012t lGrootll201 lh can al- 
low a measurement of the white dwarf spin-orbit alignment 
and rotation velocity. With these methods, this new group of 
systems will fill a poorly-covered range of the white-dwarf 
and M-dwarf mass/radius relations. 
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